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Abstract 

We perform an exact renormalization-group analysis of one-dimensional 4-state clock models with complex 
interactions. Our aim is to provide a simple explicit illustration of the behavior of the renormalization- 
group flow in a system exhibiting a rich phase diagram. In particular we study the flow in the vicinity of 
phase transitions with a first-order character, a matter that has been controversial for years. We observe 
that the flow is continuous and single- valued, even on the phase transition surface, provided that the 
renormalized Hamiltonian exist. The characteristics of such a flow are in agreement with the Nienhuis- 
Nauenberg standard scenario, and in disagreement with the "discontinuity scenario" proposed by some 
authors and recently disproved by van Enter, Fernandez and Sokal for a large class of models (with real 
interactions). However, there are some points in the space of interactions for which a renormalized Hamil- 
tonian cannot be defined. This pathological behavior is similar, and in some sense complementary, to the 
one pointed out by Griffiths, Pearce and Israel for Ising models. We explicitly see that if the transforma- 
tion is truncated so as to preserve a Hamiltonian description, the resulting flow becomes discontinuous 
and multivalued at some of these points. This suggests a possible explanation for the numerical results 
that motivated the "discontinuity scenario" . 
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1 Introduction 



One of the major features of the Renormalization Group (RG) theory is that it makes possible the 
description of the singular critical behavior associated to second-order phase transitions in terms of smooth 
RG transformations on a suitable space of Hamiltonians. By suitable we mean local (in some mathematical 
sense that must be specified) and translational-invariant Hamiltonians. The thermodynamic singularities 
associated with critical phenomena arise after an infinite number of RG steps, in the vicinity of a RG 
fixed point 0, |. 

The smoothness of the RG flow is also expected to hold for other values of the coupling constants. In 
particular, it was conjectured by Nienhuis and Nauenberg |^ that it can also be smooth around points 
where first-order phase transitions occur. In their picture [standard scenario) these points are in the 
domain of attraction of certain fixed points — called discontinuity fixed points — that govern the behavior 
of the system on the coexistence manifold. These fixed points have relevant directions whose critical 
exponents are equal to the dimensionality of the model y — D. The singularities associated to first-order 
phase transitions are obtained by an infinite iteration of RG transformations around the discontinuity 
fixed points ^, ||, in complete analogy with the case of second-order phase transitions. 

On the other hand different authors [|[ 0, ||, ^ |l^ have exhibited numerical evidence and argu- 
ments indicating that, in disagreement with the picture advocated by the standard scenario, the RG 
transformation could be discontinuous at the transition points and could associate different renormalized 
Hamiltonians to the different phases coexisting at the transition (multivaluedness). 

Recently, van Enter-Fernandez-Sokal |l^ have rigorously shown that the second picture cannot 
hold for a large class of systems. For classical variables taking values on a compact manifold with real- 
valued local interactions (e.g. absolutely summable Hamiltonians), they prove that — whenever defined — 
a renormalization-group map associated to a local (real space) RG prescription is single-valued and 
continuous on a suitable space of Hamiltonians. However, they have also shown that at — or in the 
vicinity of — a transition surface, a renormalized local Hamiltonian may not exist at all. This pathology 
is called non-Gibbsianness, and was previously pointed out by Israel [l3| and, in some sense, previously 
by Griffiths and Pearce Its occurrence has been rigorously verified |l^ for the Ising model in 
dimensions D > 2 at sufficient low temperature for some RG prescriptions (wich include decimation, 
some cases of majority rule, block-averaging and Kadanoff transformations). 

These results suggest that the discontinuities observed in Refs. (|, 0, || are in fact an artifact of the 
truncation of the renormalization scheme. If the renormalized local Hamiltonian exists, the size of the 
observed discontinuity should decrease for smaller truncations. On the other hand, it is plausible that 
the numerical manifestation of non-Gibbsianness be phase-dependent yielding an apparent discontinuity 
that persists (possibly acquiring an oscillatory character) for successively smaller truncations [ p2| . In 
relation to the validity of these explanations we mention the work of Ref . [|l5j , where truncation errors 
were estimated for Monte Carlo RG calculations of the two-dimensional Ising model below the critical 
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temperature. They were found to be of the same order of magnitude as the observed discontinuity. 

The observations of the preceding paragraph do not apply to the resuhs of Ref. JTo[ , which are not 
obtained via a truncation scheme. However, the method used there rehes on the hypothesis that there 
is only one renormalized trajectory flowing away from the discontinuity fixed point JlGf . The reported 
discontinuity of the renormalization flow could hence be a consequence of the fact that this hypothesis is 
not valid for the two-dimensional Ising model which has two relevant operators: one associated with the 
temperature {— 1/0) and the other with an external magnetic field p5[ |. 

In this paper we intend to clarify further the above issues by analyzing very simple one-dimensional 
models that on the one hand exhibit a rich phase diagram and on the other hand admit a RG scheme 
that can be exactly computed. The models belong to a family of one-dimensional g-state clock models 
with complex- valued Hamiltonians previously introduced by some of us . These models were intended 
as simplifications of quantum spin models of present interest in which cither the complex interactions 
are present ab initio or they appear in the effective Hamiltonian upon integration of fermionic degrees 
of freedom. Complex couplings arise, for instance, in the classical spin model associated to a quantmn 



Heisenberg chain ||18|, |19|, the effective model for the quantum Hall effect considered in ||2lj, 22 1, and 
the chiral Potts model which has been a recent focus of interest |2^, Q and is not unrelated to the model 
we analyze. 

On the other hand, the presence of complex interactions introduces very interesting modifications in 
the beautiful picture developed for traditional (real-interaction) statistical mechanics — pure phases, er- 
godic decomposition of states, independence of the boundary conditions for the thermodynamic potentials 
[ p5| , p6[ . This could merit further analysis, given the role played by complex interactions ||2^, |2^, |2^, ^ |l] 
in the study of several phenomena addressed by equilibrium statistical mechanics: Lee- Yang singularity 



pa, Eg, M , metastability effects Pg, P3|, high-T analyticity |P4, 35, Sq] and low-T smoothness [Bg| of 



the free energy, and deformations of the phase diagram at low temperature |33 . An important change 
incorporated by complex interactions is, of course, the enrichment of the phase diagram: the nearest- 
neighbor one-dimensional models studied here do present phase transitions, a feature totally absent in 
real-interaction models of comparable simplicity . 

The models with complex interactions that we consider here have several desirable properties |]l7| : 
They possess a Hermitean transfer matrix which is positive-definite for a wide range of coupling constants 
that includes phase transition points. This last property makes these models useful to generate unitary 
quantum- mechanical systems in the continuum limit [^8[ Although the models have complex inter- 
actions the free and internal energies are real for periodic boundary conditions. In addition, these models 
present a nontrivial phase diagram characterized by a manifold of points where the leading eigenvalues 
of the transfer matrix cross. At these transition points, the energy density is in general discontinuous, 
so the transition can be catalogued as first-order; but at the same time for many values of the couplings 
constants the correlation length diverges with a critical exponent v = I. Therefore, the surface formed 
by these points is also a critical surface. On the other hand, inside the transition manifolds there run 



curves where the thermodynamic behavior is, in some sense, even more singular. They correspond to 
the points of the phase diagram where the partition function is zero for (a sequence of) arbitrarily large 
volumes, and suitable boundary conditions. For lack of better or established nomenclature, we shall call 
these points — which have no counterpart in the phase diagrams of real interactions — "Lee- Yang- type" 
(LYT) points. This type of singularities is familiar to people studying metastability [ p2| |3^ . In principle, 
they correspond to singularities for the finite-volume free energies and its existence is not enough to in- 
fere the non-analyticity of the infinite-volume free energy^ — it only rules out the usual analyticity proof 
pionered by Lee and Yang p^ . For real interactions, the question of whether there is a singularity in 
the infinite-volume limit is related to the possibility of analytically continuing the (infinite-volume) free 
energy in the presence of metastable states [^^. In our model, however, the LYT singularities do have 
infinite-volume consequences: they belong to transition surfaces and in addition, as we shall see, the free 
energy presents some further singular behavior at these points. 

We believe that all these attributes make the clock models a useful laboratory for the study of RG flows. 
Since the models are one-dimensional and have nearest-neighbor interactions the natural candidate for 
the RG prescription is decimation. Such transformations have been used for many one-dimensional Ising 
models with real local couplings, see e.g. ]40[ |. However, these were models without first-order transition 
points. The first application to a model with this type of transitions has been presented recently The 
fiow of the decimation transformation is comparatively simple for these models: it involves a finite number 
of parameters, so it can be computed exactly without any truncation approximation. We present here 
the results for the 4-state clock model subjected to decimations with blocks of size even. The calculations 
reveal a number of instructive features which, we think, bear some light on present controversies on the 
properties of RG transformations. Let us summarize our main observations: 

1) Our results are in agreement with the standard scenario of Nienhuis and Nauenberg and the van 
Enter- Fernandez-Sokal theorems (even when, rigorously speaking, our model does not satisfy the hypoth- 
esis of these theorems because it has complex interactions): Whenever defined, the fiow is continuous and 
the behavior of the model on the surface of first-order transitions is determined by (three) discontinuity 
fixed points lying on this surface. One of these fixed points has two relevant directions characterized by 
the same critical exponent y = D = 1, as predicted by the standard scenario. The other two discontinuity 
fixed points are non-Gibbsian, and hence the fiow is singular. 

2) In some regions of the phase diagram we observe pathologies similar to those pointed out by Grif- 
fiths, Pearce and Israel |l^: Already for the first renormalization step the renormalized Hamiltonian 
fails to exist. The reason for this non-Gibbsianness is, however, different from the one observed for real 
interactions. While for the latter the renormalized measure fails to be (quasi)local ^ |l3j — roughly 
speaking it looks as if the coupling constants proliferate in an uncontrollable manner — in the present 
case the renormalized measure gives zero weight to (open) sets of configurations. This corresponds to 

^Simple example: The singular functions fn{x) = In (l — ^) converge to the perfectly analytical zero 
function. We thank Alan Sokal for clarifications regarding this point. 
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one of the coupling constants attaining the value +00. In a more abstract language, the renormalized 
measure fails to satisfy uniform non-nuUness. Such a pathology cannot happen in the renormalization of 
real interactions, except if the renormalization prescription itself excludes some configurations. Moreover, 
for real interactions quasilocality and uniform nonnuUness are necessary and sufficient conditions for a 



measure to be Gibbsian |42 (in the complex case the matter is much more involved). In this sense, we 
can say that the pathologies observed in the present example are complementary to the ones discussed 
in Refs. |l^, We observe that if at some of these points the renormalized Hamiltonian is 

"truncated" by ignoring the coupling constant that acquires an infinite value, the resulting transfor- 
mation becomes discontinuous and multivalued. In this sense we could say that in this example the 
non-Gibbsianness produces a "discontinuity scenario" as a result of truncation. 

3) We can distinguish two different "degrees" of pathological behavior. In most of the pathological 
points, the non-Gibbsianness brought by the RG transformation is "recoverable" : a further iteration of 
the RG transformation restores the Gibbsianness, and this Gibbsianness is preserved under additional 
iterations. Equivalently, the points are pathological for one particular even-block decimation scheme 
but not pathological for all the others. In constrast, at the points LYT — and only there — the non- 
Gibbsianness is not restored by a further renormalization. The points LYT are, therefore, pathological 
for all even-block decimation schemes. We observe that the LYT singularities correspond to points where 
all the eigenvalues of the transfer matrix are doubly degenerated (in absolute value). While we do not 
fully understand the meaning of this observation, we remark that a similar degeneracy seems to be present 
in the critical curve of the two-dimensional Ising model |Q . 

4) The RG flow exhibits the following characteristics: It involves two families of models, both of 
them formed by 4-state clock models parametrized by three real numbers denoted J, Ji and e, but 
differing in the value of a discrete parameter m that gives an imaginary part to the couplings. One family 
is a (transfer-matrix-Hermiteanness-preserving) complex extension of the other. For each decimation 
prescription, each family can be divided in two sections: an open region of the parameter space that 
we call the black-hole section, and its complement, the non-black-hole section. The RG transformation 
maps models in the non-black-hole section into models of the same family, and models in the black-hole 
section into models of the other family. The points where the RG transformation has a "recoverable" 
pathology are precisely those of the boundary between black-hole and non-black-hole sections that are 
not LYT singularities. These boundaries change with the decimation scheme; the only points common 
to all of them are curves of LYT singularities. 

5) Each family of models exhibits seven different fixed points (plus periodic repetitions), but three 
of them attract the majority of the models with finite couplings: a "high-temperature" (zero correlation 
length) fixed point, and the two non-Gibbsian discontinuity fixed points located on the critical surface. 
Models outside the critical surface are attracted by the high-temperature fixed point — of the same or 
the other family, depending on which side of the critical surface the initial models are. On the other 
hand, models on the critical surface are attracted by the non-Gibbsian critical fixed points of the same 



or the other family, depending on wether initiaUy the models are inside or outside the black-hole region. 
There is only a curve of points inside the critical surface that are attracted to the Nienhuis-Nauenberg 
fixed point mentioned above. Furthermore, we can explain the "terminal" pathologies affecting the points 
LYT in terms of strongly attracting invariant planes of points; By a single renormalization transformation 
the points LYT are mapped either into a non-Gibbsian critical fixed point or into an invariant plane of 
non-Gibbsian models. 

The paper is organized as follows. In Section 2 we define the 4-state clock model and study its phase 
structure in detail. The RG analysis is mainly explained in Section 3. We see that in some regions of the 
interaction space an analytic continuation of the RG equations is needed (passing to another Riemann 
sheet). This is performed in Section 4 and consists in adding an additional parameter to the Hamiltonian. 
The pathologies observed are discussed in Section 5. In Section 6 we present some final comments. 

2 Phase Structure of the 4-State Clock Model 

A one-dimensional q-state clock model is defined by a classical variable s„ ("spin") fluctuating among 
the q roots of the unity 

cos ( ) ,sin ( ) ) ; p„ = 0, . . . ,g - 1 (1) 
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on the sites n of a one-dimensional chain. We consider the Hamiltonian 

N 
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n— l,r— 1 ^ 



/ 2r7r ,\ . . f Irn , 



Jr_iCOS -(p„ -p„+i)j -f ie,. sin -(p„ -p„+i) ) \, (2) 



with Jr and Sr real parameters, and where [q/2] denotes the integer part oiq/2. This type of Hamiltonian 
is a natural generalization of the Hamiltonians studied in [ pT| , with additional interactions introduced 
so to have a system of couplings closed under decimation transformations, that is, with as many inde- 
pendent couplings (including the free energy normalization) as RG equations In the case q — ^ the 
Hamiltonian is 

N 

- I3TL = I J cos ij^jPn -Pn+l)) + Jl C0s(7r(p„ -pn+l)) 



IE sm 



-p«+l))} ■ (3) 

(We have written J, e instead of Jq^Ei to conform with previous works 0.) The transfer matrix of 



this model reads 



T . 



and its eigenvalues are 



^ie-J-i gJ+Ji ^-ie-Ji ^-J+Ji 
Q-J+Ji gie-Ji gJ+Ji g-ie-Ji 
y g-ie-Ji g-J+Ji gie-Ji g-'+'^i J 

Ao = 2e"'^ cosh J -I- 2e^^'^ cose 
A2 = 2e"'^ cosh J — 2e^-'i cose 
Ai = 2e'^^ sinh J + 2e~'^i sin e 
A3 = 2e-'^ sinh J — 2e~'-^^ sine . 



(4) 



(5) 



The eigenvectors of the transfer matrix are spin waves: the eigenvector for the eigenvalue Afc has compo- 
nents 

a;W=e"W2 (™^o,...,3). (6) 
From (||)-(||) it is easy to obtain exphcit expresions for arbitrary powers of the transfer matrix: 

3 

{T-U = (1/4) Y.^\kT exp[z^M9 - 90/2] • (7) 

With this formula we can calculate all the statistical mechanical and thermodynamic properties of the 
model. 

Before discussing the phase diagram we observe the following symmetries of the eigenvalues (||): 

Ao,2(V2-e) - A2.o(e + ^/2), Ai,3(7r/2 - e) = Ai,3(£ + 7r/2) (8) 
Ao,2(-e) = Ao,2(e), Kz{-£) = XzAe) (9) 
Ao,2(-J) = Ao.2(J), Ai,3(-J) = -A3,i(J) (10) 

Due to them, we only need and will describe the phase diagram in the region 0<e<7r/2,J>0. 
The diagram on 7r/2 < e < tt can be obtained by reflections on the plane e = 7r/2 together with the 
interchange Aq < — > A2 (symmetry (H)); the diagram on — tt < e < by reflections on the plane e — 
plus the interchange Ai < — > A3 (symmetry (^)); and the diagram for other values of e follows from 
the 27r-periodicity. The diagram for J < is obtained by reflections on the plane J = followed by 
the interchange Ai < — > —A3 (symmetry (p^)). As an example, Fig. |l] shows the dependence of those 
eigenvalues on the coupling e of the imaginary term for some choice of J and Ji . 

In the thermodynamic limit [N — > 00) the behavior of the system is regulated by the leading eigen- 
values. By "leading" we mean largest in absolute value. For those values of the parameters for which 
there is only one leading eigenvalue, say A^^, we can see from (^ that the free energy density for even 
chains takes the value^ loglA^pl and all the expectations are independent of the boundary conditions. 
By all accounts, these correspond to regions with only one phase present. Crossing points of two leading 
eigenvalues correspond to transition points, (see Fig. |l|). 

In the region < e < 7r/2, J > 0, the transition surface is composed of three pieces (Fig. |^): 

• The surface Cq.i defined by the condition Ao = Ai > max(A2, IA3I): 

Co,i = |(J, Ji,e) : < e < 7r/2 , J > ande = ^ +arcsin ^^-^=— ^1 . (11) 

• The surface Ci__3 defined by the condition Ai ~ —A3 > max (Aq, A2): 

Ci-3 = |(J, Ji,e) : 7r/4 < £ < 7r/2 , J = and £ > ^ + arcsin I . (12) 



^This also holds for odd chains with an appropiate branch choice in the logarithm definition. 
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• The surface Co, 2 defined by the condition Ao = A2 > max(|Ai|, IA3I): 

Co,2 = {(J, Ji,£) : £ = 7r/2and0 < J < 2Ji}. 



(13) 



In addition, the degenerated planes J = 00 and Ji — 00 can also be considered transition surfaces, but 
we are more interested in the properties of models with finite couplings. 

It is trivial to see that the energy density is discontinuous at the transition points, which implies that 
the system undergoes a first-order phase transition. However, the correlation length also diverges at the 
points of the surface Co,i and this transition surface can also be considered as a critical surface. This 
can be seen, for instance, by computing the correlation of two spin variables 



where Xko is the leading eigenvalue. Eq. ( |l4| ) implies that the correlation length diverges at the transition 
points Co,i as ^(/3) — {f3^^ — f3^^)^^ and hence one obtains critical exponents i' — rj — 1. The correlation 
length also diverges when J goes to infinity and, in fact, in the limit J — > cx) the discontinuity of the 
energy density goes to zero and the transition becomes purely second order. 

An important feature brought by the use of complex couplings is the presence of "Lee- Yang-type" 
(LYT) singularities. We so denominate those values of the couplings for which there exists a sequence of 
volumes and boundary conditions giving a partition function equal to zero. For the present models this 
condition requires the existence of a sequence of powers of the transfer matrix with some entry (which 
could be different for different powers) equal to zero. Therefore, the locus of the points LYT can easily 
be obtained from (0) together with the form (||) of the eigenvalues A^. The conclusion is that, for finite 
couplings, there are only two curves of points LYT (again, for the region < e < 7r/2, J > 0; for the rest 
of the phase diagram one must proceed as commented below (|l0|)): 





(14) 




• The line defined by the equations Aq = A2 and Ai — —A3; that is: 



LYTl = {(J, Ji,e) : J = 0,£ = tt/2}. 



(15) 



For couplings on this line, 




n is even and \q — g'| = 1, 3 . 



(16) 



• The curve defined by the equations Aq = Ai and A2 = —A3; that is: 



LYT2 = {(J, Ji,e) e Co.i : cosh2J cos2£ 



1} 



(17) 



On this curve. 




n is even and |g — (/'j — 2 



(18) 
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(See Fig. We remark that these are the only solutions if we search for a sequence of volumes with 
zero partition functions. In addition there are whole surfaces in the coupling-constant space formed by 
points for which the partition function is zero for some volume and boundary condition. These points 
are precisely those of the boundary dB2k of the black-hole region (see next section) for some even-block 
decimation. We note, in passing, that the diagonal entries {T'^"')qq are always nonvanishing, in fact 
strictly positive (for finite couplings). 

Both curves LYT are contained in transition surfaces: LYT2 is in Co,i and LYTl is in Co, 2 for Ji > 
and in Ci,_3 for Ji < 0. In fact, the curves correspond to those points of the transition manifold where 
all the eigenvalues become pairwise degenerate (in absolute value) . The curves intersect at the point 



which correspond to all the eigenvalues having the same absolute value (maximun degeneration). This 
point is extremely singular for the flow of the decimation transformations (Section 

Within the critical manifold Co,i we observe some thermodynamic features very different from those 
of real interactions: the infinite-volume free energy acquires a dependence on the boundary conditions, 
and no boundary condition produces truncated correlations with decaying behaviour. Let us denote fq,q' 
the free energy for boundary conditions q on the left and q' on the right. From (^ we obtain, for example, 
that in the surface Co,i, fq^q' — logAo (=logAi) for all boundary conditions q,q' except if q' ~ q ± 2. 
In that case the free energy takes the values logA2 or logjAsj depending on which is the subleading 
eigenvalue. The transition between these two values takes place at the curve LYT2 — both non-leading 
eigenvalues become equal in absolute value — where fq^q±2 can only be defined by considering volumes of 
odd size. 

For the construction of quantum-mechanical systems via scaling limit, only the region of parameters 
for which the model satisfies the reflection-positivity condition is of interest. For odd chains, this condition 
is satisfied for all values of the parameters because the transfer matrix is Hermitean . For even chains 
the reflection-positivity holds only if all the eigenvalues of the transfer matrix are non-negative, that is, 
for J > J]ip{ Ji,s), with 



(Figs. ^ and ^). In this work we do not restrict ourselves to the reflection-positive region; if we did 
we would miss the curves LYT which are the most interesting regions from the point of view of the 
pathologies of the decimation transformation. 



Plyt = ( J = , Ji = , e = tt/2) 



(19) 




(20) 
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3 Renormalization Group Flow 



A single one-spin decimation transformation in the system described above yields the following Renormalization- 
Group equations 

1, ri + e4"'icosh2J 
J = o log 



2 ^ [ e^'-'i + cos 2e 

, _ 1 / (1 + e-"^ cos2£)(l + e^"^^ cosh2J) \ 

"^1 " 4^°^\ 2(cosli2J + cos2£) ^ 

. 1 f 1 + cosh2Jcos2e 

e — - arccos 



2 [ cosh 2 J + cos 2e 

where the primes denote the renormalized quantities. To decide the quadrant of e' , the last equation 
must be complemented with 

tane' = tanetanhJ. (22) 

This is the complete exact RG flow except that, as usual, we have omitted the renormalization of the 
coupling associated with the identity operator. This corresponds to a spin-independent constant added to 
the Hamiltonian, which is therefore irrelevant for the analysis of correlation functions and thermodynamic 
potentials. By iteration of ( pT[ ) we can understand the flow of all decimation transformations of even 
blocks (decimation of an odd number of spins = even powers of T). 

Starting with the coarser features of the flow, we first observe that all such transformations map the 
region J < Jrp{Ji, s) (non-reflection-positive region) onto the complementary region J > Jrp{Ji, s) by 
a single RG transformation [Q. This is due to the fact that the matrix is always reflection positive. 
In particular all the flxed points lie in the reflection-positive region. We then notice that the flow defined 
by (|l|)-(|22|) exhibits a periodicity of tt in e and very distinct symmetries with respect to the change of 
sign of J and e. We conclude that it is enough to study the flow in the region of Fig. ^ — ( J > and 
< e < ■"'/S) — and so we will in the sequel. The flow for J > and — 7r/2 < e < is an identical copy, 
while the points with J < change quadrant: their renormalized couplings are the same as those for | J|, 
but with opposite sign for e' . Having said this, we can forget about the quadrant-fixing equation ([2^). 

Another conspicuous feature of the RG transformation (^ij) is the presence of a large region in the 
coupling constant space where the renormalized couplings are not real numbers. We shall call this region 
the black hole B2 (see Figs. | and |): 



So 



|(J, Ji,e) : e*"'! < -cos2e, e > ^1 (23) 



Notice that the intersection of the black hole region with the transition surface is non-empty, even in the 
region where reflection positivity holds (Fig. |3|). While ( p3| ) corresponds to the decimation of alternated 
spins, decimations of chains of k spins also exhibit a corresponding black- hole region B2k- In Section 
^ we discuss how this region changes with the order k of the decimation. In the black-hole region we 
must extend analytically ( ^l| ) into the complex space (J, Ji,£) G C^. However, this extension must be 
performed so as to preserve a physically highly desirable property that motivated the initial restriction 



to real (J, Ji,e): a real free energy. For this we extend ( |2l| ) in such a way that the transfer matrix T 
remains Hermitean. This extension will be done in detail in Section ^ 
The boundary 



of the black hole is a pathological region for the RG flow: at these points the renormalized couplings 
take the value J' — +00, J{ = —00. As we discuss in more detail in Section this implies that the 
renormalized measure gives zero weight to some sets of configurations, a property incompatible with a 
Boltzman probability weight (exponentials are never zero). The renormalized Hamiltonian simply does 
not exist. There is another pathological region, namely the line LYTl ( J — 0, e = tt/2), where J( = +00 
and e' is undefined. In Section || we shall discuss more extensively these pathologies. 

In the rest of the coupling constant space (including the points on the transition surface not in dB2 or 
LYTl), the flow is well defined: to every point the transformation (pl| ) associates a unique renormalized 
Hamiltonian. Such a non-pathological flow will be described in the remaining part of this section. 

Due to the special features of the decimation procedure we can easily locate RG invariant surfaces, 
that is, surfaces that are mapped onto themselves by a RG transformation. One of them is the reflection 
positivity interface defined by (|2^) which is defined by the condition that the smallest eigenvalue of the 
transfer matrix be equal to zero. This condition is invariant under a decimation transformation, a fact 
that can be more easily seen if we write this transformation in the form 



where L is the (strictly positive) overall factor associated to the renormalization of the free energy. The 
same expression (|25| ) shows that equalities and inequalities among the absolute values of eigenvalues of 
T are also preserved by decimations of any even block, and, therefore, that the refiection-positive part of 
the transition surfaces are invariant under these transformations. In particular, the surfaces Co,i, Co, 2, 
J = 00 and Ji = 00 are RG invariant. The intersections of two invariant surfaces are RG invariant 
curves, that is, renormalized trajectories. Some of them are drawn in Fig. |^a. 

These renormalized trajectories intersect at the RG fixed points (see Figs. || and ^). They are listed 
m Table |l|. From (|l|) it can be seen that the correlation length diverges at Fc, F*, i^g^ and Fq. 




(24) 



XliJ, Ji, e) = Afe( J', J[,e') L{J', J[,e'), 



(25) 
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Point Position (J, Ji, e) 



Type 



F 


(0,0,0) 


stable 


Ft 


(0,00,0) 

(i logs, i logs, ^/2) 


saddle point 


Fc 


saddle point 


F* 


(oo, — oo, 7r/4) 


saddle point 


Ft 

^ C 


(0,00, V2) 


saddle point 


F"t 


(oo, oo, e) ; e e [0, n/4]) 


unstable 


^^0^ 


(oo, ilogcos2£,£) ; e g [0, 7r/4]) 


saddle point 



Table 1: Fixed points for even-block decimations. 



The linearization of the RG equations (21) around the different fixed points yields the matrices 




11 i tan 2e 
1 i -itan2£ 







Lpet 



-| 
2 




(26) 



1 



whose eigenvalues (Q) and eigenvectors determine the critical exponents (y, = log Ci/ log 2) and the main 
properties of the renormalized trajectories. 

The fixed point Fqo is completely attractive in all directions. The point F^ behaves as an attractor 
with respect to perturbations along the J-axis and the e-axis, but has a marginal direction (y = 0) along 
the Ji-axis. The latter is actually a relevant direction, as can be checked numerically, and flows toward 
the infinite-temperature fixed point Fqq. 

There are three fixed points on the transition manifold: Fc, FJ and F*. The first one has features 
predicted by the standard scenario: It is a saddle point with one irrelevant and two relevant directions, 
both with critical exponents equal to y = D = I. Of the relevant directions, one corresponds to perturba- 
tions along the e-axis with the flow escaping towards the flxed point Fqo , while the other is tangent to the 
renormalized trajectory J — Jb,p{Ji, ""/S) where the flow moves towards Fj and F*. The flow towards F* 
is along the renormalized trajectory defined by the intersection between the reflection-positivity interface 
and the transition surface Co,i (given by ( |20| ) and ( pT| ) respectively); the flow towards Ft is along the 
line J = Jj^p( Ji, 7r/2). The irrelevant direction is along the RG trajectory J — 2Ji, e — tt/2. 

The other two flxed points on the transition surface correspond to singularities of the flow (they 
are non-Gibbsian) and therefore the RG transformations can not be linearized around them. We have 



11 



analyzed the flow numerically and we conclude that is attractive for those transition points satisfying 
J < 2Ji (triangular- looking cusp in the surface of Fig. ||a), while F* attracts the rest of points of the 
transition surface. In both fixed points there is one relevant direction that flows away from the transition 
surface, whose critical exponent is not well defined. 

Finally we have two ("zero-temperature") fixed-point lines Fq\ and Fq. There are three marginal 
directions at each point of Fq^, but only the one tangent to the line Fq^ itself is truly marginal. The 
other two are actually relevant and the fixed points become unstable. Each of the points of Fq has 
one irrelevant perturbation {y = —1) along the (2,-1,0) direction and two marginal ones: one, truly 
marginal, tangent to the line Fq itself, and the other given by the lines Ji = oo, e e [0, 7r/4], which are 
RG trajectories fiowing from Fq'' to Fq. At these fixed points one can take the continuum limit. The 
existence of these lines of fixed points leads to different quantum mechanical systems in the scaling limit 




We now analyze the flow outside of the black- hole section B2, to follow the points inside the black 
hole we need an analytic extension to be discussed in Section ^ (See Figs. ^ and ||b). For reading the 
rest of this section we recommend to look at the Figures ^, |^a, ^a and ^ as many times as needed. We 
start with the flow inside the invariant planes. There are four of them, corresponding to e = 0, e = 7r/2, 
Ji — 00 and J — 00. From ( |2l] ) we see that F^o attracts all points in the invariant plane e = 0, except 
those sitting at the line Ji = 00 which flow to F^. The fixed points Fq^ are completely unstable. In the 
plane e — t:/2 the situation is less simple: Fc attracts the transition points on the renormalized trajectory 
J — 2Ji and FJ is attractive for all transition points above this line (J < 2Ji). The remaining points in 
this plane flow toward the black hole region in a finite number of RG steps, except those of the boundary 
dB2 which go to the plane J — 00 (at its intersection with Ji — —00) in one renormalization step. In 
the plane Ji = 00 all points are attracted by Fj^, except those belonging to the line e — tt/2 which flow 
toward Fj. We remark that the points lying on the line J — are mapped to F^ in a single RG step. 
Finally, in the plane J = 00 every line e = constant is a trajectory of the renormalization transformation. 
When e G [7r/4,7r/2] they flow towards the black hole region, while the points with e G [0, 7r/4) are 
attracted by Fq and repeled by Fq^ 

In the rest of the space the flow is as follows. The points above the critical surface Cq.i (see Fig. |^) 
are driven by the RG transformation to the infinite-temperature fixed point F^o. They will reach this 
point after an infinite number of RG steps. In particular, the points at the plane J — are mapped into 
the plane e = in a single renormalization step, where they remain during their migration to Foe- (The 
points with J < go to the previous quadrant — 7r/2 < e < 0.) The points below Co.i are swallowed up 
by the black hole region in a finite number of RG steps. That means that the part of the black-hole region 
B2k located below the surface Co,i grows with the order k of the decimation procedure. For instance, the 
corresponding black-hole region for the transformation T ^ (double decimation) contains, besides 
all the points belonging to B2, also the points which are mapped into the black hole B2 by a single 
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decimation. (In Section ^ we shall see that the part of i?2fc above the surface Co,i decreases in size with 
k.) Finally, the points on the transition surface Cq.i remain there after renormalization. There are three 
domains of attraction: a) if J = 2Ji (e = tt/2) the flow goes towards Fc; b) if J < 2Ji (e ~ 7r/2), towards 
F^; and c) if J > 2Ji, toward F^. In all cases, an infinite number of RG steps are needed to arrive at 
the corresponding fixed point. 

4 Analytic Extensions of the 4-State Clock Model 



As explained in the last section, at every point belonging to the region B2 one has to extend (21) 
analytically to the complex {J,Ji,e) space^. However, to have a real free energy this extension must 
be performed in such a way that the transfer matrix remains Hermitean. This amounts to a suitable 
definition of the complex logarithm. Alternatively, we can ask ourselves how many hermitean extensions 
of a matrix T of the form (||) can be constructed. After some straightforward algebra we conclude that 
there are two of them, which can be labelled by an integer m taking the values and 1. The transfer 
matrix now reads 



T(to) 



/ g./+./i g-ie-Ji ^_^^m^-J+Ji gie-Ji \ 



gie— Ji gJ+Ji p-ie — Ji 



pie-./i g./+./i g-i£-Ji 



(-l)™e"''+-'i e 



(27) 



where {J,Ji,e) £ R^. This implies that the couplings which appear in the Hamiltonian (|^) must be 
replaced by the following quantities 



J ^ j_ H|ZL 



T T 1 trmr 

^ + 4 (28) 



'1 'Ji 

e ^ e 



A ^ A+i^ 

where A is the real coupling associated to the identity operator. Notice that e always remains real. The 
case m — corresponds to the original clock model (^), and m = 1 to its analytic continuation with 
coupling constants carrying an imaginary part. This last model will be called the extended clock model 
and quantities related with it will be denoted with a bar. In Ref. pl[ | the parameter m was not considered. 

In order to study the phase structure of this extended model, one has to compute the eigenvalues of 
the new transfer matrix T = r(l) [given in (p7|)]: 

Ao = 2e''i sinh J + 2e~'^^ cos e 

A2 = 26"^^ sinh J — 2e^''^ cose , . 

Ai = 2e''i cosh J + 2e"'^i sine ^ ' 

A3 = 2e'^i cosh J + 2e"-'i sine. 

We can relate easily these ones with those given by (H) 

A™(J, Ji,e) = A™+i(J, Ji,e + 7r/2) . (30) 
^ We thank Alexei Morozov for pointing out the interest of this question. 
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This relation implies that the phase structure of the extended model can be obtained from that of the 
original one (m — 0) by reflecting Fig. ^a with respect to the plane e = -Kj^. Thus, the RG flow in the 
extended model is the same as in the original one, except for a reflection (Fig. ||b). In particular, to each 
fixed point located at (J, Ji, e) in the original model (listed in Section 3), corresponds an analogous fixed 
point but located at (J, Ji, 7r/2 — e). For the same reason there exists another black hole 

B2 = {(J, Ji,e) : e*'^i < cos 2£, e < ^ ; m = 1} (31) 

where an analytic extension of the corresponding RG equations is needed. The curves of LYT singularities 
are now 



and 



which intersect at 



LYTl = {(J, Ji,e) : J = = e} (32) 
I:yT2 = {{J,Ji,e) e Co,i : cosh2J cos2e = 1} , (33) 
Plyt = (J = 0,Ji = 0,£ = 0). (34) 



The new feature is the flow of points inside the black-hole regions B2 and B2 from one model to the 
other. To study it we need the RG equations in the slightly higher-dimensional space (J, Ji, s, m), for a 
transfer matrix T{m) of the form (p7|). These are 

J' = llog((-l)-+-' l + ^"^-^»^h2J 



2 ' e4-fi + (-l)™cos2e 

_ 1 / (1 + (-l)"'e-4^^ cos2£)(l e^-^^ cosh2J) ) 

" 4^°n^ ^ 2(cosh2J+(-l)-cos2£) / ^^^^ 

, 1 r (-1)'" +cosh2Jcos2e1 

£ = 7: arccos < — — ; r — > , 

2 \cosh2J+ (-l)'"cos2£ J 

where the renormalized integer m! must be chosen in such a way that there exist real numbers J' and J( 
satisfying (|35|). This choice is unique. 

Consider the original clock model (m = 0). If (J, Ji,e) does not belong to the black hole B2, then it 
follows from (^5|) that ml = 0. The flow in this region was described in detail in Section 3. However if 
(J, Ji,s) G B2 then one must choose m' = 1; so the renormalized Hamiltonian belongs to the extended 
model. By the same arguments, if to = 1 the renormalized Hamiltonian will have to' = 1 whenever we 
start at a point not belonging to B2, but m' = if (J, Ji,s) G B2. Thus, the analytic extension of the 
extended model is the original one. 

Furthermore, the decimation transformation sends the points inside a black hole to a different side of 
the critical surface: A point (J, Ji, e, to = 0) G ^2 placed below Co,i is mapped to a point (J', J[, e', to' — 
1) ^ B2 located above Co,i. Such a point is eventually attracted by the fixed point -Fqo, except if it 
belongs to the line J = 00 in which case is attracted by the line Fg (see Figs. ^ and ^). On the other 
hand, if (J, Ji,e, m = 0) belongs to the part of the black hole on top of Co,i, it is mapped to a point 
outside B2 and below Co,i- Such a point ends up, after a finite number of iterations, inside B2 and below 

14 



Co,i, so it is returned to the m = model outside the black hole B2 and above Co,i. The attractor for 
such a point is therefore Foo- 

In addition, a point on the critical surface Co,i and inside B2 goes to a point on the critical surface 
Co,i of the extended model. This image point is outside B2, so it is eventually attracted by F* 

In particular, all this implies that when the order k of the decimation increases, the part of the black 
hole i?2fc on top of Co.i decreases in size, with its boundary approaching Cq.i. Complementarily, the 
part of i?2fc located below Cq.i increases in size, with the boundary also approaching Cq^i for large k. 
Therefore, in the limit of decimations of large order k, the black-hole region occupies the whole region 
situated below Co,i, and the boundary dB2k asymptotically coincides with this critical surface. All these 



boundaries 3i?2fc contain the curve LYT2 defined in (17), and they "turn" around it when they approach 
Co,i as fc ^ 00. Moreover, the curve LYT2 is the intersection of each 9i?2fc with Co,i and it divides 9i?2fc 
in two parts: a part inside the black hole B2k (above the curve LYT2 in Fig. ^), and a part outside it. 
In Section |^ we discuss other important features of the RG transformation at this curve. 

Analogous considerations hold for the extended (m ~ 1) model. In particular, points in B2 and below 
Cq.i are attracted by the fixed points Fqo or Fq of the original model; while points in B2 and above Co,i 
are attracted by F^o- The points in B2 n Co,i are attracted by F*. 

Finally, the points on the transition surface Ci._3 are attracted by the fixed point f]., except those 
belonging to the curve Ci,_3 n Co,i which fiow to Fc- 

This concludes the analysis of the flow at the non-pathological points, both in the original and the 
extended model. We observe no trace of ambiguities: to each Hamiltonian there corresponds a unique 
renormalized Hamiltonian given by (p5[). Besides, all fixed points are reached after an infinite number of 
RG steps, except F^^ and which act as very strong attractors (attract in one RG step) for the points 
on the lines J = 0, Ji = cxd, to = and J = 0, Ji = 00, m = 1 respectively. 

What is left is the study of the pathological regions: the boundaries dB2 and dB2 — or, more generally, 
dB2k and dB2k — of the black-hole regions, and the line LYTl. This is the subject of next section. 

5 Pathologies in the 4-State Clock Model 

Griffiths and Pearce |lj] were, to our knowledge, the first to point out possible "peculiarities" for RG 
transformations at a first-order phase transition. This pioneer call of attention was followed by numerical 
results and tentative arguments ^ ||, H, |l^ that seemed to indicate the possibility of discontinuity 
and multivaluedness of the flow at the coexistence points. The theorems proved in ]l^ have ruled out 
such lack of smoothness for real interactions and compact single-spin space, and have left the more subtle 



phenomenon of non-Gibbsianness |11, ^ ^ as the main pathology for real-space renormalization. 

For real interactions, the class of Gibbsian measures is exactly characterized by the class of quasilocal 
and uniformly non-null measures . Roughly speaking, quasilocality means that the direct influence of 
far away spins on a given spin s^^ is small. That is, if we take a set A centered in xq and fix the spins 



15 



Sx G A, a; 7^ a;o, then the change of the expected value of with the boundary conditions outside A 
is vanishingly small as the diameter of A goes to infinity. On the other hand, a necessary condition for 
a measure to be uniformly non-null is that every open set of configurations have nonzero measure. For 
example, the measures obtained for spin systems in the zero-temperature limit are typically not uniformly 
non-null as they are concentrated on the configurations which minimize the energy of the system. For 
general complex interactions the theory of Gibbs measures is not well developed, but in any case the 
exponential form of the Boltzman weights implies that uniform non-nuUness is also a necessary condition 
for Gibbsianness. 

The non-Gibbsianness of the renormalized measure exhibited in ^ |l^ for Ising models at or close 
to a first-order phase transition, is a consequence of lack of quasilocality. In the present case, however, the 
origin of the pathologies is rather related to the loss of uniform non-nuUness. The pathologies occur when 
there are renormalized models for which some matrix elements of the transfer matrix (^7|) vanishes. For 
instance, if the transfer matrix is diagonal it would imply that the corresponding measure is concentrated 
on constant configurations. In terms of coupling constants, measures that are not uniformly-non-null 
correspond to manifolds where J or Ji take an infinite value, and the pathological points are those 
mapped into such manifolds by a finite number of iterations of the RG transformation. 

We remark that the "pathological" character of a RG transformation at these points only appears if 
we insist on finding a renormalized Hamiltonian. The renormalized transfer matrix is always well defined. 
This remark is equivalent to the fact, emphasized in ||ll|, that renormalization transformations (for 
real interactions) are always well-defined as maps between probability measures. It is only at the level of 
Hamiltonians that the induced transformation can become sick. 

To see at which points the single-spin decimation is pathological, we have to look at the matrix 
elements of obtained from (|2^). These are: 

iT^)q,q = 2 (e^-'i cosh2J + e-2.^i) 

(r2),,,+i = 2(e-^-- + (-l)"e-'^+-) = (r2);,+3 (36) 

iT^)q,q+2 - 2((-l)™e2'^i+e-2^icos2£). 

Zeros of any of these elements determine points that are pathological for a single RG map. The diagonal 
elements are always non- vanishing; however, the off-diagonal ones can vanish at some points of the 
coupling constant space. On the one hand, {T'^)q,q+2 = at the boundaries dB2 and dB2 of the black 
hole region, and {T^)q,q+i = {T^)q.q+3 = at the lines LYTl and LYTl. The matrix T'^ is diagonal at 
Plyt and Plyt- 

It is interesting to follow the flow of the points at these pathological regions. For concreteness we 
analyze the original (m — 0) model; the discussion for the extended model is analogous. We see that 
all the points in 9i?2 sue mapped into the line J' = oo, J{ = — oo: The points above LYT2 go onto 
the segment < e' < 7r/4 (i.e. outside the black-hole region B2), while those below LYT2 end up at 
the complementary segment 7r/4 < e' < tt/2 (deep inside ^2)- AH the points in LYT2 are mapped to 
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the point with e' = 7r/4, i.e. to the non-Gibbsian fixed point F*. Applying a further iteration of the 
transformation one can see (using, for instance, the expresion for obtained via (|^), that the points 
with e' 7^ 7r/4 recover finitely- valued couplings: Those with < e' < 7r/4 are mapped to points above 
Co.i and flow towards Foo; while those with 7r/4 < e' < 7r/2 are mapped to points above Co,i (m = 1 
model) and flow towards F^c- On the other hand, the points in LYT2 remain in F* (it is a fixed point !) 
and so they lead always to non-Gibbsian renormalized measures. This non-Gibbsianness is the result of 
the strong attraction that a non-Gibbsian fixed point has on the curve LYT2 (in the sense of attracting 
it in finitely many — in fact one — steps). 

The points LYTl are also victims of a strong attraction: In one decimation transformation the 
half-line J = 0, Ji > 0,e = n/2 is mapped into the plane Ji = oo, while the other half-line — 
J — 0, Ji < 0, e = tt/2 — is mapped into the plane Ji = oo. As these planes are RG-invariant and 
formed by non-Gibbsian models, those points remain non-Gibbsian under iterations of rcnormalization 
group transformations, ft is suggestive to observe that, at the level of coupling constants, the RG trans- 
formations exhibit discontinuities and multivaluedness at LYTl. For example, points arbitrarily close 
to the line LYTl and contained in the plane J = are renormalized into points with e' = 0, while 
equally close points but contained in the plane e — t:/2 remain with e' = tt/2. This shows an explicit 



discontinuity. Moreover, the value of e' given in (21) is undefined at the line LYTl: By approaching this 
line ("preparing the system") from different directions J(e) one can obtain any value of e'. This is an 
instance of multivaluedness. In this sense, the "most" pathological point is Plyt = LYTl n LYT2. Its 
renormalized transfer matrix is diagonal and, moreover, both e' and J{ are undefined for this point, and 
can be made to take any value by approaching it in different ways. 

Again we emphasize that this discontinuity and multivaluedness appears because we are trying to 
follow finitely-valued couplings, and "truncating" the infinitely-valued ones. There is no ambiguity or 
lack of smoothness at the level of renormalized transfer matrices. Analogously, non-Gibbsianness prevents 
us from linearizing the transformation as a function of the couplings around the critical fixed points F* 
and F^ , and hence the matrices Lp^ and Lpt are not defined. The linearization is perfectly possible for 
the transformation as a function of the entries of the transfer matrix. 

Other pathologies of the single-spin decimation show up if wc look for zeroes among the entries of 
higher powers of T. Such pathologies correspond to the loss of Gibbsianness after a finite number of 
renormalizations. They coincide with single-renormalization pathologies for decimations of higher order, 
and they occur at the successive boundaries 9i?2fe- Such pathologies have the same features as those 
discussed above: except at the points LYT2 — contained in all such boundaries — the Gibbsianness is 
recovered for good after a further iteration of the transformation {dB2k H dB2k' = LYT2 \i k ^ k'). On 
the other hand the above pathological points depend very much on the rcnormalization group prescription 
and in this sense thery are not universal. However their existence in the espace of local Hamiltonians 
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cannot be avoided by the choice of a suitable choice of a renormaHzation group prescription]^ In fact 
the analysis in a general decimation prescription can be obtained from the study of the continuous flow 
defined by the infinitesimal transformation T ^ T* {t E R). This is given (for m = 0) by the following 
differential equations 

T / X , ^ » , ^ N cosh J , , , ^ , , ^ , sinh J 
J = (Ai logAi + A3 logAsj^^j^^ (AologAo + A2logA2) ^^j^ 

■ _ (Aq logAp + A2 logA2) (Aq logAp - A2 logA2) . tane - tanhJ 
^ " Se'^coshJ Se-'^icose ^2 2 ^ ' 

£ = — {(Ai logAi - A3 log A3) cose - (AologAo - A2 log A2) sine} 

which are valid only in the reflection-positive region. These equations are well defined for every point 
in this domain which contains some points of (?i?2- The only pathology of the points of B2 is that they 
are attracted by the line J — 00, Ji — —00, e £ (7r/4, 7r/2] in a finite "time" t < 2. This feature explains 
the especial behavior of the points of B2 under the action of T^-renormalization group transformations: 
they are mapped into the to = 1 domain. The singularity arises also in this infinitesimal scheme when we 
integrate ( p7| ) up to a finite renormalization scale t. For each t there is a black hole region Bt that changes 
as we increase the scale t. In fact, as we have seen, every point below the transition surface belongs to a 
region Bt for t large enough because they are attacted by the line J — 00, Ji — —00, e G (vr/4, tt/2). This 
implies that the flow defined by the differential equations ( |37| ) does not generate a one-parameter group 
of global transformations. 

6 Conclusions 

We believe that the model studied here provides an instructive illustration of the behavior of real-space 
renormalization transformations in the vicinity of coexistence manifolds. Admittedly, some of the features 
observed could be specifically due to the use in particular of decimation procedures, but we expect others 
to be indications of phenomena of more general validity. 

The present example explicitly exhibits a picture that takes considerable more work to verify for 
real interactions |pT| , [l^ , p^ : the reconciliation of the standard scenario and the presence of pathologies, 
and non-Gibbsianness as a source of an observed discontinuity and multivaluedness in the flow of the 
coupling constants. In addition, we consider singularly suggestive the observations that the worst type of 
pathologies takes place at points where all the eigenvalues of the transfer matrix are doubly degenerate, 
and where there are singularities of the type linked, in some studies, to metastability effects. 

The analysis performed here can be extended to g-state clock models of higher q. We do not expect 
any new phenomenon except, perhaps, in the limit q 00 (chain of plane rotors with complex interac- 
tions). But of course, the challenge is to find a model with real interactions exhibiting, with comparable 



^In some cases a description of the renormalization group transformation in terms of non-local variables 
might lead to a continous flow Q . 
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explicitness, non-Gibbsianness of renormalizcd measures and its consequences for computational schemes 
(an aspect left unresolved in previous work 13, Such an example, however, may not be easy to 
construct given the rather subtle character of non-quasilocality. 



Acknowledgements 

We thank Antonio Gonzalez- Arroyo and Alexei Morozov for useful discussions, Aernout C. D. van Enter 
for carefully reading and critizicing the manuscript, and Prof. Michael Fisher for correspondence. We 
also thank to Prof. Joel Lebowitz for bringing Ref. ||3l[ to our attention. R. Fernandez and J. Salas thank 
the members of Departamento de Fisica Teorica of the Universidad de Zaragoza and of the Universidad 
Autonoma de Madrid (R.F.) for their warm hospitality during the completion of this work. We also 
acknowledge CICyT for partial finantial support (grants AEN90-0029 and AEN90-0030). 



19 



References 

[1] K. G. Wilson, Rev. Mod. Phys. 47 (1975) 773. 

[2] K. G. Wilson and J. Kogut, Phys. Repts. 12 (1974) 75. 

[3] B. Nienhuis and M. Nauenberg, Phys. Rev. Lett. 35, (1975) 477. 

[4] W. Klein, D. J. Wallace and R. K. P. Zia, Phys. Rev. Lett. 37, (1976) 639. 

[5] M. E. Fisher and A. N. Berker, Phys. Rev. B26, (1982) 2507. 

[6] H. J. Blote and R. H. Swendsen, Phys. Rev. Lett. 43 (1979) 799. 

[7] C. B. Lang, Nucl. Phys. B280 [FS18] (1987) 255. 

[8] A. Gonzalez- Arroyo, M. Okawa and Y. Shimizu, Phys. Rev. Lett. 60 (1988) 487. 

[9] A. Hasenfratz and R Hasenfratz, Nucl. Phys. B295 [FS21] (1988) 1. 
[10] K. Decker, A. Hasenfratz and P. Hasenfratz, Nucl. Phys. B295 [FS21] (1988) 21. 
[11] A. C. D. van Enter, R. Fernandez and A. Sokal, Phys. Rev. Lett. 66 (1991) 3253. 
[12] A. C. D. van Enter, R. Fernandez and A. Sokal, in preparation. 

[13] R. B. Israel, in Random Fields Vol. II, Eds. J. Fritz, J. L. Lebowitz, D. Szasz. North Holland, 
Amsterdam (1981). 

[14] R. B. Griffiths and P. A. Pearce, Phys. Rev. Lett. 41 (1978) 917; J. Stat. Phys. 20 (1979) 449. 
[15] A. Gonzalez-Arroyo and J. Salas, Phys. Lett. B261 (1991) 415. 

[16] K. G. Wilson, in Recent Developments in Gauge Theories, eds. G.'t Hooft et al. Plenum, New York 
(1980). 

[17] M. Asorey and J. G. Esteve, J. Stat. Phys. 65 (1991) 483. 

[18] I. Affleck, Nucl. Phys. B257 [FS14] (1985) 397. 

[19] F. D. Haldane, Phys. Rev. Lett. 50 (1983) 1153. 

[20] I. Affleck, Phys. Rev. B43 (1991) 3215. 

[21] A. M. M. Pruisken, Phys. Rev. B15 (1985) 2636. 

[22] H. Levine and S. B. Libby, Phy. Lett. B150 (1985) 182. 

20 



[23] H. Au-Yang, B. M. McCoy, J. H. H. Perk, S. Tang and M. L. Young, Phys. Lett. A123 (1987) 219. 
[24] R. J. Baxter, J. H. H. Perk and H. Au-Yang, Phys. Lett. A128 (1988) 138. 

[25] D. Ruelle, Statistical Mechanics, Benjamin (1969); Thermodynamic Formalism, Addison Wesley 
(1978). 

[26] R. B. Israel, Convexity in the Theory of Lattice Gases, Princeton University Press (1979). 

[27] T. D. Lee and C. N. Yang, Phys. Rev. 87 (1952) 404; 410. 

[28] P. J. Kortman and R. B. Griffiths, Phys. Rev. Lett. 27 (1971) 1439. 

[29] M. E. Fisher, Phys. Rev. Lett. 40 (1978) 1610; Prog. Theor. Phys. Suppl. 69 (1980) 14. 

[30] J. L. Cardy, Phys. Rev. Lett. 54 (1985) 1354. 

[31] G. Gallavotti, J. Lebowitz, Physica 70 (1973) 219. 

[32] S. N. Isakov, Comm. Math. Phys. 95 (1984) 443. 

[33] M. Zahradnik, J. Stat. Phys. 47 (1987) 725. 

[34] G. Gallavotti and S. Miracle-Sole, Comm. Math. Phys. 7 (1968) 274. 
[35] R. B. Israel, Comm. Math. Phys. 50 (1976) 245. 

[36] R. L. Dobrushin and M. R. Martirosyan, Teor. Mat. Fiz. 74 (1988) 16. 

[37] H. O. Georgii, Gihhs Measures and Phase Transitions, de Gruyter, Berlin-New York (1988) and 
references therein. 

[38] M. Asorey, J.G. Esteve and A. F. Pacheco, Phys. Rev D27 (1983) 1852. 
[39] C. Borgs and J. Z. Imbrie, Comm. Math. Phys. 123 (1989) 305. 
[40] M. E. Fisher, D. R. Nelson, Ann. Phys. (N. Y.) 91 (1975) 226. 
[41] M. Asorey, J. G. Esteve and J. Salas, DFTUZ 91.32 preprint. 
[42] O. K. Kozlov, Prohl. Inform. Transmission 10 (1974) 258. 

[43] R. A. Minlos and Ya. G. Sinai, Theor. Mat. Fiz. 2 (1970) 230; translation Theor. Math. Phys. 2 
(1971) 167. 

[44] K. Gawedzki, R. Kotecky and A. Kupiainen, Comm. Math. Phys. 47 (1987) 701. 

21 



Figure Captions 

Figure 1: Eigenvalues (fc = 0, . . . , 3) of the transfer matrix (5) as a funetions of e for J — 1.5 and 
Ji = 0.5. 

Figure 2: Transition surfaces and LYT points for the model in the region J > 0,0 < e < 7r/2. The 
intersection with the reflection-positivity interface of Fig. ^ is also shown. 

Figure 3: Reflection positivity interface (|2^) for the model. We have also drawn some of the fixed points, 
the renormalized trajectory joining Fc and F* and the black hole region (^). 

Figure 4: Transition surface outside the black hole region (^). The contour of the reflection positivity 
interface is also drawn for clarity. 

Figure 5: Rcnormalization group flow in the reflection positive region: (a) for the original model, and 
(b) for the extended one. Fixed points and some of the renormalized trajectories are also despicted. The 
black hole region is despicted in the right side of (a) and in the left one of (b). 

Figure 6: Universality classes (a) for e G (0,7r/2), and (b) e — 7r/2. The fixed point which attracts all 
points belonging to a region is written inside it. Notice that the surface Co, -3, defined by the condition 
Aq = —A3 > max(|Ai|, IA2I), is a critical one in the region J < (See Eq. (10)). The dashed line 
represents the boundary dB2 and the dash-dotted line the LYTl points. 
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